Transition from weak to strong cascade in MHD turbulence 
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The transition from weak to strong turbulence when passing from large to small scales in magnetohydrody- 
namic (MHD) turbulence with guide field is a cornerstone of anisotropic turbulence theory. We present the first 
check of this transition, using the Shell-RMHD which combines a shell model of perpendicular nonlinear cou- 
pling and linear propagation along the guide field. This model allows us to reach Reynolds numbers around 10 6 . 
We obtain surprisingly good agreement with the theoretical predictions, with a reduced perpendicular energy 
spectrum scaling as k'] 2 at large scales and as fc~ 5/3 at small scales, where critical balance between nonlinear and 
propagation time is reached. However, even in the strong regime, a high level of excitation is found in the weak 
coupling region of Fourier space, which is due to the rich frequency spectrum of large eddies. A corollary is that 
the reduced parallel spectral slope is not a definite test of the spectral anisotropy, contrary to standard belief. 
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Introduction. — Plasma turbulence plays an important role 
in solar corona, solar wind, fusion devices and interstellar 
medium. The knowledge of turbulent energy spectra is a basic 
step to solving problems like cosmic ray transport, the tur- 
bulent dynamo and solar corona or solar wind heating. In 
this letter, we focus on incompressible magnetohydrodynam- 
ics (MHD) turbulence with strong guide field. A first theory 
has been first proposed in the limit of small relative magnetic 
fluctuation b/Bo, assuming isotropy of the cascade HHH]. This 
theory is based on the weak coupling of Alfven waves, the 
Alfven decorrelation time t a decreasing more rapidly than the 
nonlinear time t^L with wave number. However, numerical 
and experimental evidence point to an anisotropic cascade, 
mainly in the directions perpendicular to the guide field fl3HZ|]. 
This prompted (0] to include the anisotropy within the defini- 
tion of the turbulence strength as x — t a /tNL — kj_b / (k\\Bo) 
where the Alfven and nonlinear times now involve respec- 
tively the parallel and perpendicular components of wave vec- 
tors. They thus proposed that the weak anisotropic cascade 
proceeds from large scales with ^ <s 1 toward smaller per- 
pendicular scales with the parallel scales remaining fixed, so 
that^f would increase and reach unity at some scale. The cas- 
cade would then be weak at large scales and strong at smaller 
scales where time scales would remain equal x — 1 (critical 
balance or CB). This change of regime should be character- 
ized by a clear-cut change of spectral scaling with wave num- 
ber k, from k~ 2 to k~ 5/3 . 

The status of the weak and strong cascades are not identi- 
cal: while an analytical approach of the weak cascade is possi- 
ble 19l[l0t], the strong cascade theory remains a phenomenol- 
ogy. The ansatz on weak/strong transition is a fundamental 
assumption in theories of anisotropic MHD turbulence lUUl . 
but it seems to contradict numerical results of Reduced MHD 
(RMHD) simulations. In RMHD nonlinear couplings along 
the guide field are suppressed, which is a valid approximation 
when the guide field is strong enough 11121 1 1 311 . When xo was 
decreased starting from 1, the spectral slope was found to vary 
smoothly, between -5/3 and -3 (with volume forcing, IU4I0 . 



or between -3/2 and -2 (with boundary forcing, [ 15111611 ). In 
no case was the spectrum found to exhibit the predicted break 
between two inertial ranges. The -3/2 spectral slope (instead 
of the -5/3 value) was ascribed to a systematic weakening of 
interactions due to a local dominance of one Elsasser species 
over the other 01711 , These results indicate that the concepts 
of strong and weak cascades might not be as robust as gener- 
ally believed, although the low Reynolds numbers reached in 
direct numerical simulations of RMHD might be the cause of 
the discrepancy. 

We thus revisit this issue here, using the shell model for 
RMHD HEl], in which the perpendicular coupling terms 
are simplified compared to RMHD, thus allowing us to reach 
Reynolds numbers Re 10 6 . In this model, the duality be- 
tween the perpendicular spectrum and the parallel space is re- 
tained, allowing us to test the transition from weak to strong 
cascade. We consider here Shell-RMHD with volume forcing 
and periodic boundaries, varying the parameter xo in order to 
compare with iflill . We will first obtain the coexistence of 
the strong (-5/3) or weak (-2) slopes along the inertial range 
when forcing large eddies with^o < 1, in complete agreement 
with CB theory. Accordingly, the eddy correlation time will 
be found to agree either with the nonlinear time (in the strong 
coupling range) or with the Alfven propagation time (in the 
weak coupling range). We will also show that high frequency 
fluctuations of the large eddies induce an unexpected excita- 
tion level in the high k\\ region of the spectrum. 

Equations. — The RMHD equations with mean field Z?o par- 
allel to Ox axis read for the variables z ± — u + b, with V.z* = 



ot ox " po 



(1) 



where viscosity and resistivity (v) are assumed equal. The 
shell-RMHD is obtained from RMHD by replacing the per- 
pendicular nonlinear couplings at each point x by a dynamical 
system defined in Fourier space. The perpendicular Fourier 
plane is paved with N + I concentric shells, the model retain- 
ing one wave number k n per shell, and one complex scalar 
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mode z„(x, t) per shell, with |z,tl 2 /2 being the total energy 
the shell n. The equations read 018111911 : 



in 



k ±n = 2"k n = O...N 
d t zt±B Q d x £ = T±-vkle n +f* 



(2) 
(3) 



where / n + = f~ = f„(x,t) represent the large scale (kinetic) 
forcing, non zero for n = 0, 1, 2 (see below), and the are the 
nonlinear terms. These are made of a discrete sum of terms of 
the form: 7** = Ak m z^Zg with m, p, q being close to n, which 
replace the convolution terms resulting from the Fourier trans- 
form with respect to the perpendicular coordinates of the orig- 
inal RMHD equations Eq.Q] The control parameter is the time 
ratio xo imposed by the forcing term f*: 



Xo 



t°/t° 

a' NL 



k[b rms /(k{.B ) 



(4) 



where k^ and k± are the characteristic wave numbers of the 
forcing respectively in the parallel and perpendicular direc- 
tions (note that b rms - u rms ). 

Numerical method and parameters. — We choose as in 
Perez and Boldyrev [14] the forcing correlation time t{ or to 
be smaller than the large scale Alfven time fi a = (k^Bo)' 1 , 
Bo = 5, and the aspect ratio of the domain to be L ± /L z = 1/5. 
We define Kq = 2n/L z and ko = 2n/L±. The smallest per- 
pendicular forced wave number is always kf ± - ky (for a spec- 
trum to develop, the Shell model needs to be excited at least 
in three consecutive shells ko, 2ko, and 4ko). Adding larger 
scales above forced scales doesn't modify the results. We 
consider two forcings, varying k^ : a narrow one or strong 

forcing (k^ = 2Kq, xo = 1), and a wide one or weak forc- 
ing (£jj e [2, 64]Kq, xo - 1/32), estimating a priori in both 

cases b rms = 1, kr x ^ 2, and = max(£jj) in Eq.|4j 

Starting from the solutions z%(x, t) = z ± (x,k ± ,t) of Eq. [3] 
(hereafter we drop the index n in k ±n ), we define three time 
scales that will be used to verify if turbulence satisfies the 
critical balance (CB) condition. We focus on the z + signal, 
and check that using z~ leads to the same results. The corre- 
lation time t cor (k ± ) is defined as the full-width-half-maximum 
(FWHM) of the autocorrelation function in time A(z + ), of the 
signal z + (x, k±, t), computed at each position x and then aver- 
aged over the spatial domain (which is homogeneous). The 
correlation length l} cor {k±) is defined in a similar way, as the 
FWHM of the autocorrelation function in space A(z + ) x com- 
puted at each time t and then averaged on 25?jvl (statistically 
stationary time series). Definitions are summarized below, as 
well as the turbulence strength x(k±Y 



tcorikj = <FWHM[A(z + ),]>, 

t a (k±) = <FWHM[A(z + ). t ]),/B„ 

tffdk±) = l/(k ± z ms (k ± )) 

X(k±) = t a (k ± )/t NL (k_L) 



Ll r {k ± )/B 



(5) 
(6) 
(7) 
(8) 



where {z) s = 1/5 J zds stands for the average, A(z) s — 
1/5 J z(s)z(s-s')ds' for the autocorrelation function, and s = 



t, x indicate the coordinate over which the average/correlation 
function are computed (time and space respectively). With the 
above definitions, the CB condition reads x = 1 or t a = t^L- 
Note that the quantity t a is the lifetime of a signal of coher- 
ence size l} cor propagating at the Alfven speed, which should 
also be equal to t cor : this will be checked as well. 

Finally we give the expressions for the energy spectra. De- 
noting by z ± (A:||, k ± , t) the Fourier transforms with respect to x 
of the signals z ± , we compute the 3D time-averaged spectra, 
£3 , the ID reduced parallel and perpendicular spectra y as 

Ei(k h k x ) = k-^mk^tt), (9) 

El(k x ) = 2n J" dkk ± E*(k h k x ) = k- 1 {\z ± (x,k ± ,t)\\ t (10) 
Efiky) = 2n j dk ± k ± E^(k h k ± ) (11) 

Note that E* = f E±d 3 k = f E±dk ± = J 

For further comparison, we define the theoretical 3D angu- 
lar spectrum resulting from the strong cascade as [HI 



E CB (h,k ± ) = k- m - q - l f(x) 



(12) 



where m is the slope of the ID perpendicular spectrum and q 
follows from the CB condition (assuming m = 5/3 we have 
q = 2/3). The boundary of the excited spectrum is given by 
the x - 1 curve, its parallel extent is given by the CB con- 
dition: k\\ oc k 2 / 3 . This is expressed by f(\x\ > 1) — 1 and 
f(\x\ < 1) - 0. Applying the definitions Eqs. (11011111) one 
obtains the familiar reduced spectra E± oc k 5 / 3 and Ef °c k^ 2 . 
Eq. [12] requires adopting a frame attached to the local propa- 
gation axis, i.e., the local mean field. However here the prop- 
agation axis is always Bq, in the absence of nonlocal terms in 
the shell model. This allows to consider eq.[T2]to be valid in 
the absolute frame. 

Results. — In Fig. Q] we show the reduced perpendicular 

total energy spectrum for the strong (xq = 1) and weak 

5/3 

Otfo = 1/64) forcing cases, compensated by fe ± . With strong 
forcing (left panel), the scaling is close to £~ 5 ' 3 in the inter- 
val 100 < k± < 10 4 . With weak forcing (right panel) one 
sees two power laws, the scaling being k~ 2 at large scales and 
Ai~ 5/3 at small scales. Long integration times (25 fi NL ) as well 
as large Reynolds number are needed to reveal this composite 
spectrum. We found equipartition between magnetic and ki- 
netic energies in the weak k~ 2 range, while magnetic energy 
dominates by a uniform factor =* 2 in the strong range (not 
shown). 

We now examine the time scales defined in Eqs. (07]). They 
are plotted versus k±_ in Fig. [2] Except at the forcing scales 
(gray band) and a bit below, the autocorrelation and the Alfven 
times (resp. dashed and solid lines) are superposed, showing 
coherence of the method. In the strong forcing case (Fig.|2^), 
the correlation time is about twice the nonlinear time (dotted 
line) in the whole inertial -5/3 range (CB condition). In the 
weak forcing case (Fig. Wp), the correlation time is constant 
at large scales, being given by the Alfven time based on the 
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FIG. 1: (a): xo = 1, (b): xo = 1/32. Reduced perpendicular (to- 
tal) energy spectra compensated by fc~ 5/3 (solid lines). Perpendicular 
forcing scales are indicated as shaded areas in each panel. Dotted 
lines: reduced spectra built from the £ 3 spectrum with^ < 1/2 exci- 
tation suppressed. The dot-dashed line is the k^ 2 scaling. 
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FIG. 2: (a): xo = 1- (b): xo — 1/32. Characteristic times t A , t cor , and 
t NL for the signal z + (x, k ± ,t) as defined in Eqs. 15171 Forcing scales 
are marked in gray. 



parallel forcing scale. Then, at small scales, it switches to a 
value which is again about twice the nonlinear time. This is 
compatible with the CB theory, which predicts the change in 
perpendicular spectral slope observed previously in Fig. [TJ), 
the CB condition holding for k ± > 10 3 . 

In Fig. [3]we show the reduced parallel total energy spec- 
tra, again for the strong and weak forcing cases (panel (a) and 
(b) respectively) compensated by k^. With strong forcing, 
the parallel spectrum scales as Hence it is steeper 

than the perpendicular spectrum but flatter than the CB pre- 
diction (slope -2). With weak forcing the spectrum at low 
wave numbers is shaped by the forced parallel modes (gray 
band on top). Their signature persists because the cascade is 
weak there (t A <sc t NL in Fig. Wp)- For higher wave numbers 
one expects to find the strong cascade scaling as in panel (a), 
on the contrary the slope is steeper, even steeper than the -2 
predicted by the CB phenomenology. 

The different parallel slopes found in the reduced spectra 
can be understood by considering the underlying 3D angu- 
lar spectrum. We first examine the strong forcing case, Fig.|4] 
Isocontours of time-averaged 3D spectrum £3 (fey, k ± ) are plot- 
ted in the left panel. We give two representations of the 
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FIG. 3: (a): xo = 1- (b): Xo = 1/32. Same format as Fig.[TJbut for 
the reduced parallel energy spectra. The dot-dashed line is the k^ 2 
scaling. 



CB: a theoretical x — 1 contour, adopting strictly the scal- 
ing k\\ oc k 2 ^ (dotted line) and an effective x — 1 contour, 
using the computed b rms (kj_) (dashed line). They differ sig- 
nificantly only when entering in the perpendicular dissipation 
range. The CB condition, ;f = 1, separates two regions: strong 
coupling is found on the left of the x — 1 l me ^ weak cou- 
pling on the right. One can see that the energy isocontours 
are mostly horizontal in the strong region ^ > 1 ; the spacing 
between successive contours corresponds nearly to the scaling 
k~ W/3 (Eq.[T3, in the range 10 2 < k ± < 10 4 . A remarkable 
feature is that, contrary to the usual assumption, the energy 
density doesn't drop abruptly when entering the weak cou- 
pling region (x < 1), at least for large perpendicular scales. 
This is best seen in the horizontal cuts (right panel of Fig. |4|, 
where again the CB is overplotted. For £3 > 10~ 13 (corre- 
sponding to k± < 10 4 ), the energy density £3 decreases as 
k^ 2 . The extent of this scaling decreases progressively with 
increasing k ± . In the CB phenomenology, the reduced parallel 
scaling in k^ 2 is due to the contribution of parallel wave num- 
bers satisfying^ = 1, which occurs only at small perpendicu- 
lar scales since parallel excitation is assumed to be negligible 
for^ < 1. Here, on the contrary, the dominant contribution to 
the reduced parallel scaling is due to the large perpendicular 
scales (with;^ <k 1). This is demonstrated in Fig. [3^ where the 
dotted line, with the typical scaling k^ 2 , represents the equiva- 
lent CB spectrum, i.e. the reduced parallel spectrum obtained 
from £3 after suppressing excitation in the region^ < 1/2: 
one sees that the resulting spectrum is much closer to the k^ 2 
scaling than the full solution. 

We consider now the weak forcing case in Fig. [5] The en- 
ergy isoncontours for k ± < 10 3 reveal a cascade at constant 
parallel wave number (k\\ < 10 2 ). Then, at higher k ± , the 
parallel spectrum widens according to the scaling k\\ oc k 2 / 3 . 
Horizontal cuts of the 3D spectrum (right panel) have features 
similar to the strong forcing case: an inertial strong cascade 
and a steep spectrum at intermediate and large perpendicular 
scales respectively, but now shifted to the last parallel forced 
mode (ifcjf a 10 2 ). A gain there is no abrupt energy decrease 
when passing the^- = 1 boundary, however the scaling is now 
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FIG. 4: Anisotropy properties (xq = 1). Left panel: contours of 
energy spectrum E 3 in the (k^,k ± ) plane. Right panel: horizontal 
cuts of £3 vs k\\ for k ± = ko 2" with n = to 16 from top to bottom. 
Dot-dashed line: k^ 2 scaling. In both panels the dotted line is the 
theoretical boundary x = 1 ; me dashed line is the measured boundary 
X = 1 • Forced scales are marked by a shaded area. 



be due to arbitrariness in defining the nonlinear time. It is 
worthwhile noting that reduced spectra in the solar wind agree 
with the form dictated by the CB condition, but that an excess 
of parallel energy may be required 112311 . consistent with our 
findings. Much is to be gained by testing turbulence theories 
using the Shell-RMHD model, especially with respect to the 
two points just mentioned: the conditions of appearance for 
the local scaling (local z + '/z~ imbalance) which controls the 
global spectral scaling and the origin of the frequency spec- 
trum of large eddies which controls the anisotropy of turbu- 
lence. 

We thank W.-C. Miiller, M. Velli, and 6. Giircan for use- 
ful discussions. A.V. acknowledges support from the Belgian 
Federal Science Policy Office (ESA-PRODEX program). 
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FIG. 5: Anisotropy properties (xo = 1/32). Same caption as Fig. [4] 
The dot-dashed line is the k^ 3 scaling. 



oc fc~ 3 , thus explaining the steeper spectral slope found in the 
reduced parallel spectrum in Fig. |3j). 

Discussion. — In Reduced MHD, parallel structures reflect 
directly the temporal structure of the perpendicular nonlinear 
excitations via the linear propagation of Alfven waves. Hence, 
the £|j~ 2 scaling of the large perpendicular scales should result 
from a f~ 2 spectrum of the low perpendicular wave numbers 
(in the weak case, the slope -3 plays the same role). We 
have checked that the f~ 2 spectrum is indeed present (see 

TM1). 
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This scaling is not due to forcing since 



suppressing nonlinear couplings leads to a steeper spectrum 
(slope A:^ 4 ). The k^ 2 scaling develops in the unforced case af- 
ter a few nonlinear times and lasts for more than ten nonlinear 
times. We conjecture that the high frequencies in large ed- 
dies are due to an inverse transfer from small to large scales, 
which may occur even in absence of a proper inverse cascade 
associated with a definite invariant I22I1 . 

Critical balance tests gave good results, with t cor about 
twice t^L in the strong inertial range (Fig. |2}. We have found 
that, although globally balanced, the two Alfven species spec- 
tral amplitudes show also a local imbalance of a factor 3. This 
might explain the factor two in time scales, which as well may 
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